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Abstract 

The kinetics of hydrogen absorption by magnesium bulk is affected by two main activated pro- 
cesses: the dissociation of the H2 molecule and the diffusion of atomic H into the bulk. In order to 
have fast absorption kinetics both activated processed need to have a low barrier. Here we report 
a systematic ab-initio density functional theory investigation of H2 dissociation and subsequent 
atomic H diffusion on TM(=Ti,V,Zr,Fe,Ru,Co,Rh,Ni,Pd,Cu,Ag)-doped Mg(OOOl) surfaces. The 
calculations show that doping the surface with TM's on the left of the periodic table eliminates 
the barrier for the dissociation of the molecule, but the H atoms bind very strongly to the TM, 
therefore hindering diffusion. Conversely, TM's on the right of the periodic table don't bind H, 
however, they do not reduce the barrier to dissociate H2 significantly. Our results show that Fe, 
Ni and Rh, and to some extent Co and Pd, are all exceptions, combining low activation barriers 
for both processes, with Ni being the best possible choice. 
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I. INTRODUCTION 



Hydrogen is regarded by many as a possible energy vehicle (or fuel) for future mobile 
applications, and has been targeted to replace the current use of liquid hydrocarbons in 
the next few decades Unlike fossil fuels, it is an environmentally friendly, non-polluting 
fuel simply because its combustion product is water, providing it is produced by renewable 
energy resources, obviously. 

One of the main challenges faced by the development of the so called hydrogen economy 
is the capability of storing hydrogen safely and efficiently. For mobile applications the 
storage materials need to satisfy a number of requirements: i) they need to be capable of 
storing hydrogen in excess of 6.5% in weight, ii) the kinetics of absorption has to be fast, 
i.e. with time scales of minutes, and iii) the temperature at which they release hydrogen 
(the decomposition temperature) needs to be ideally in the range 20-100°C. Cyclability of 
the material is also a desired property. 

Metal hydrides are natural hydrogen storing materials, and the relatively strong bonds 
between hydrogen and the host metal satisfy the safety requirement. Unfortunately, however, 
no material which has all the properties mentioned above exists today (see for example the 
review by Sakintuna et al. [2I). 

Magnesium hydride, MgH2, satisfies some of the above requirements. It has high storage 
capacity (7.6 wt %), good cyclability and it is relatively inexpensive. However, its enthalpy of 
formation is too high (-76 kJ/mol), requiring temperatures in excess of 300°C to decompose 
it into H2 and Mg bulk. The formation of the hydride also has slow kinetics js, 4|, making 
this material not good enough. However, MgH2 represents a good test material to study how 
various treatments can affect its properties. In particular, it has been found that by doping 
the material with transition metals can weaken the Mg-H bond and reduce the stability 
of the MHde (see fo. example fl fl fl fl fl U. and .efe.nces tKe.in). Ban-.U„„g 
further enhance the sorption processes by increasing the number of possible paths for the 

diffusion of H (see for example fly 

. \\\ and references therein). A new method of chemical 
:luid deposition in supercritical fluids (SCF) has also recently been used on metal hydrides 
1^ . This method offers the same sorption properties of ball milled samples, but with a 
hugely improved cyclability (the catalytic effect of the metal being almost constant for SCF 
samples, while decreasing for ball-milled samples after about 100 cycles) and in particular 
shows that the catalytic effect of Ni on hydrogen sorption processes is higher than that of 
Pd. 



Liang et al. [l^ and Schulz, Liang & Huot [l^ found that V and Ti are better catalysts 
than Ni for hydrogen absorption and desorption from MgH2-metal composites, showing faster 
absorption kinetics at T ~ 300 °C and faster desorption kinetics above ~ 250 °C than other 
3(i-elements investigated. They also found an enthalpy of hydride formation for different 
catalysts similar to that of MgH^. By contrast, theoretical calculations and experimental 

n n 

results of Song et al. [7| and Shang et al. [8| led to different conclusions. They found that 
the stability of MgH2-Ni is reduced when compared to that of MgH2-Ti. Moreover (and 
contrary to the experimental findings of Liang et al. jisl), the heat of formation of the 
metal-doped MgH2 hydrides is smaller than that of MgH2. In particular, MgNiH2 shows a 
smaller enthalpy of formation than MgTiH2. 

Zaluska et al. [sl used Li, Al, V, Mn, Zr and Y as catalysts for the hydrogena- 
tion/dehydrogenation of Mg alloy samples. According to their results, V remarkably im- 
proves the H absorption kinetics, but Zr is better for lower temperature H. However, the 
best kinetic results are achieved with mixtures, i.e. V + Zr or Mn + Zr Mg alloys. Bo- 
bet et a.. Q have sKown tKat the M-gen storage ptope.t,es are enhanced when using 
reactive mechanical alloying of Mg + 10 wt.% Co, Ni and Fe mixtures. Co, unlike Ni, is 
found to significantly increase the quantity of MgH2 formed. However, Bobet et al. later 
reported that the hydrogen sorption properties of Mg-Co mixtures are less effective than 
those reported for MgH2-metal mixtures. Gutfleish et al. [l^ have recently presented re- 
sults achieved with a Mg sample alloyed with Ni (1 wt %) and Pd (0.2 wt %). Their sample 
shows excellent hydrogen absorption/desorption kinetics and cyclic stability, exhibiting an 
overall reversible H2 storage capacity of 6.3 wt. %. 

Previous theoretical and experimental investigations over pure surfaces of transition met- 
als belonging to the left of the periodic table have shown that H2 dissociation is promoted, 
but also that the bonding between the hydrogen atoms and the metal is strong (see [l7| and 
references therein; see also isl. Iiol). 

In our previous paper 20|] we have shown that H2 dissociation on the metal(Ni,Ti)-doped 
Mg surface has a barrier similar to that on the corresponding pure metal (111) surface, 
and also that the strength of the hydrogen-metal bond is similar to that on the pure metal 
surface. The strength of the H-TM (TM=Ni,Ti) bonding was found to be correlated to the 
height of the diffusion barrier. We therefore might expect to see an analogous trend in the 
dissociation and diffusion barriers by doping the Mg(OOOl) surface with various transition 
metals. In fact, we will show that the elements on the left of the periodic table make the 



H2 dissociation barrier to vanish but are responsible for high diffusion barriers, while those 
on the right cannot catalyse the dissociation of the molecule. Among the elements studied 
here, we found that Fe, Ni and Rh, and to some extent Co and Pd offer a good compromise 
between the promotion of dissociation and the hindering of diffusion, and qualify as good 
catalysts for accelerating the kinetics of hydrogen absorption. 



II. COMPUTATIONAL METHOD 



All the DFT calculations were performed with the ab-initio simulation package VASP 



using the projector augmented wave (PAW) method 22, l23| and the PBE exchange- 



2l| 



correlation functional 
up the calculations 



2J]. An efficient charge density extrapolation was used to speed 



. We used a plane-wave basis set to expand the electronic wave- 
functions with the same plane- wave energy cut-off of 270 eV as in [2^, which guarantees 
convergence of adsorption energies within 1 meV. Surfaces were modeled using periodic slabs, 
with 5 atomic layers and a vacuum thickness of 10 A. The topmost three atomic layers were 
allowed to relax, while the bottom two were held fixed to the positions of bulk Mg. Calcu- 
lations were performed using 2x2 surface unit cells, with 9x9x1 k-point grids and replacing 
one of the four surface Mg atoms by one TM(=Ti,V,Zr,Fe,Ru,Co,Rh,Ni,Pd,Cu,Ag) atom. 
These settings were extensively tested and guarantee convergence of activation energies to 
better than 0.02 eV. Activation energies have been calculated with the nudged elastic band 



(NEB) method [26|] using 17 replicas, which proved to be sufficient to reach convergence of 
activation energies to better than 0.01 eV, and display all the main features of the minimum 
energy path. The initial state of the NEB calculations for the dissociation of H2 is repre- 
sented by the hydrogen molecule sitting on top the TM at a distance of 5 A, and the final 
state is the most energetically favourable among four possible adsorption sites for the two 
dissociated hydrogen atoms (see Ref. [2^ for details). For the diffusion process, the initial 
state is represented by the final state of the dissociation process, and the final state by a 
configuration where one of the two hydrogen atoms has been displaced into a nearby hollow 
site (see details in Section [III Bl) . 



Figs. [2] and [3] have been made using the XCRYSDEN software 27 1. 



III. RESULTS 



In the following section we report calculations for the bulk structural parameters of the 
various elements investigated, as a test of the quality of the PAW and the PBE exchange- 
correlation functionals. In Sec. IIIIBI we report results for the H2 dissociation and diffusion 
barriers, which we also analyse in terms of the position of the centre of the d-band of the 
various transition metals employed as dopants on the Mg(OOOl) surface. 



A. Bulk parameters 

We obtained bulk structural properties of the pure transition metals by calculating en- 



ergy versus volume curves, and fitting them to a Birch-Murnaghan equation of state 28 1 



The elements investigated here were: Zr, V, Fe, Ru, Co, Rh, Pd, Cu and Ag, together with 



Ti and Ni already presented in Ref. 20|]. The bulk parameters were derived using 13x13x13 
and 18x18x12 k-point grids for those metals with the cubic and the hexagonal structure 
respectively. The corresponding standard version of the PAW functional was used for all of 
them, with the exception of V and Zr for which we used the version of the PAW treating 
respectively the 3s^3p^3d^4s^ and 4s^4p^4d^5s^ electrons in valence (which give results closer 
to the experimental values than the standard versions of the PAW functionals). Calculated 
bulk parameters values are reported in Table [U together with the details of the PAW poten- 
tials, which include the electrons treated in valence and the core radii. Overall, the lattice 
parameter a is always overestimated, and the bulk modulus is underestimated with respect 
to the experimental values. This is in agreement with the findings from previous theoretical 
calculations, and is typical of the PBE functional. Among the elements investigated, only 



Ni, Co and Fe are magnetic. We find a magnetic moment of 0.63 20|, 1.70 and 2.15 /x^/atom 



for Ni, Co and Fe respectively, which are in agreement with the corresponding experimental 



values of 0.61, 1.71 and 2.22 //^/atom 3l|. However, as discussed in the next Section, Fe is 



the only element which required spin-polarised calculations. 



B. H2 dissociation and diffusion 



The activation barriers for H2 dissociation over the various metal-doped Mg surfaces are 
reported in Table [Tll where we also report two experimental values for the H2 dissocia- 
tion/recombination on the Mg(OOOl) surface. The value reported in Ref. 48|] (~ 1.0 eV) 



TABLE I: Bulk properties of pure transition metals (TM=Ti,V,Fe,Co,Ni,Cu,Zr,Ru,Rh,Pd,Ag). For 
each element we report the bulk lattice constant a (together with c/a for hep metals), the bulk 
modulus kg, the electrons treated as valence (VE) and the core radius Tcore of the PAW potentials. 



References for values previously reported in literature follow in the last column. 
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TABLE II: Activation energies for II2 dissociation (E^iss) on the pure Mg and metal-doped Mg 
surfaces (ordered by increasing atomic number). 



Metal surface Ediss (eV) 

pure Mg 0.87", 0.4^'^ 0.5'^'^ 1.15^, 1.059, 0.95 

[Expt.] 1.0*, 0.75±0.15J' 
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from DFT LDA calculations and PES. This lower value as compared to other calculations is 
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refers to the recombination barrier (which, in this particular case, is similar to the dissocia- 
tion barrier) identified with the barrier for desorption from the surface. This value was not 
directly measured in the thermal programmed desorption (TPD) experiments of Ref. 48 1 
because complete desorption spectra as function of temperature could not be taken, due to 



the onset of Mg sublimation at ~ 450 K which overlaps with the temperature at which H2 
desorbs. However, it was noted that the onset of H2 desorption appears at 425 K, which is 
similar to that of the H/Be(0001) system that has a determined desorption energy of ~ 1 
eV {5^, and so, by analogy, it was suggested that the activation energy for desorption might 
be the same on the H/Mg(0001) system too. The value reported in Ref. 49|] (0.75±0.15 eV) 
has been obtained by the interpretation of TPD experiments performed on a 400 A thick 
magnesium film. This value is in good agreement with the calculated PBE dissociation 
energy, but is significantly lower than the dissociation energies calculated with RPBE by 



Vegge 



45] and Du et al. 46|] which are 1.15 and 1.05 eV respectively. It should be noted. 



however, as pointed out in Ref. 49|], that the experimental situation may not be the same 
as the theoretical ones, due to the possible presence of steps on the surface which might be 
more reactive sites and lower the H2 dissociation barrier. Moreover, the inferred dissociation 
energy of 0.75±0.15 eV is based on the use of the Arrhenius relations with assumed pre- 
factors of ~ 10^^ Hz. As showed in Refs. 5JJ] and 52], these values could be underestimated 
by more than two orders of magnitudes because the classical pre-factors do not include the 
enhancement due to the large entropy increase as the molecules leave the surface, in which 
case the activation energy could be up to ~ 0.25 eV higher. 

The geometry of adsorption of the H atoms on the metal-doped Mg surfaces listed in 
Table [III appears to be somewhat correlated to the height of the dissociation barrier. We 
find that when the barrier is large (i.e., for Cu, Pd, Ag) the H atoms fall into filled hollow 
sites, and when the barrier is null (i.e., for Ti, Zr, V and Ru) they fall into empty hollow 
sites. In between there are elements showing a small energy barrier for which the preference 
towards filled hollow sites (i.e., for Ni) instead of empty hollow sites (i.e., for Fe, Co and 
Rh) is weaker (less than about 30 meV). 

Fe is the only dopant for which magnetic calculations are really required, with the total 
magnetic moment on the Fe atom being of 2.8 and 2.5 fiB in the initial and final state of 
the dissociation process respectively (Co is magnetic too, but when used as dopant of the 
Mg surface our calculations show that it can be treated as non-magnetic). In particular, 
non-magnetic calculations for the Fe-doped surface would give significantly different results, 
reducing the dissociation barrier to almost zero and increasing by 60 % the energy difference 
between the initial and final states. 

As we noted before 20|, the activation barrier for the dissociation of H2 over a metal- 



doped Mg surface is similar to that on the corresponding pure metal surface. For example. 



in the case of Cu our calculated barrier is 0.56 eV, which is close to the dissociation barrier 
of a 



53 



Dout 0.5 eV suggested by experiments and other DFT (GGA and PBE) calculations [19|, 



54 



55 



56 



57( 1 over the pure Cu(lll) surface. In addition, we find that the energy 



difference between the final state and the initial state is -0.19 eV, which is in line with the 

who found that H2 dissociation on the pure Cu(lll) surface 



findings of Kratzer et al. 



is exothermic, with a gain of 0.2 eV. On the Pd-doped surface we calculate a dissociation 
barrier of 0.39 eV, which reduces to 0.30 eV when a smaller 5x5x1 k-points grid is used 
(instead of the 9x9x1 grid). The reason for mentioning the result obtained with the coarser 
grid is because we want to compare with the findings of Dong, Kresse & Hafner 59|], who 
performed calculations with a similar grid for the dissociation of H2 on the Pd(lll) surface, 
and found a barrier of 0.29 eV when the molecule dissociates on top a Pd atom, therefore 
a value very close to our value of 0.30 eV on the Pd-doped Mg surface (note, however, 
that they report the bridge-bridge as the preferred dissociation path, with a barrier of only 
70 meV, which is consistent with the theoretical value found by Nobuhara et al. 6^ and 
in good agreement with the experimental value of 50 meV 6l|)- Finally, in the case of 
Ag, our activation barrier of 1.18 eV is in agreement with the experimental results which 
predict a dissociation barrier on Ag(lll) larger than that on Cu and somewhat larger than 
0.8 eV and we also agree with previous theoretical calculations 63] also obtained on 
the pure Ag(lll) surface, which reported an activation energy of 1.11 eV. Fig. [T] shows the 
MEP's for hydrogen dissociation on a TM= V, Fe, Co, Cu, Zr, Ru, Rh, Pd, Ag doped Mg 
surface investigated here, together with those on a TM=Ni, Ti -doped Mg surface that we 
have reported previously [20 1. 

Figure [2] shows the dissociation of H2 on the Ag-doped Mg surface, which is found to 
have the largest barrier value among all the dopants investigated here (since this barrier is 
larger than that on pure Mg, obviously H2 will not dissociate onto this site, but will rather 
choose some regions of the Mg surface free of Ag). For Fe, Co, Cu, Rh and Pd doped Mg 
surfaces, the images at the IS, TS and FS of the NEB are similar, with the H2 molecule at 
the TS sitting closer (i.e., Cu and Pd) or further away from the surface (i.e., Fe, Co and 
Rh) compared to the behavior shown in our previous paper 20|] for the Ni-doped Mg surface 
(see dn-surf values reported in Table [In|) . For Ag we note that the hydrogen molecule at 
the TS dissociates closer to the surface than on the Ni-doped Mg surface, and that it does 
so on a side of the dopant atom (see Figure [2]). 

A closer look at the geometry of the dissociation process shows an interesting correlation 



TABLE III: The d-band center position with respect to the Fermi energy (Ed), the activation 
energy barrier for the dissociation of H2 (E'diss)^ the energy difference between the final and initial 
state i^di^^^^) dissociation, the activation energy barrier for the diffusion of atomic H 

(Ediff) and the corresponding energy difference CE^ffJ^^) on the pure Mg surface as opposed to 
the metal-doped Mg surfaces (these have been ordered so as to highlight the overall dependence 
along each column of the periodic table, as we go from right to left across the periodic table). 
Also reported in the last column is the average distance of molecular hydrogen from the surface as 
measured at the transition state, dH_surf- 
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between the height of the barrier and the geometry of the transition state. The dissociation 
of molecular hydrogen into two hydrogen atoms happens on-top of the dopant atom when Fe, 
Co, Ni or Rh are used as dopants, slightly shifted to the side when Pd and Cu are the metals 
used as dopants, and fully on the side when Ag is the dopant. In other words, it appears 
that H2 dissociates on top of the dopant atom for those doped Mg surfaces which show a 
very small barrier (i.e., Fe, Co, Ni and Rh), slightly shifted to the side of the dopant atom 
for the Pd and Cu doped Mg surfaces having a non negligible barrier, and completely on the 
side of the dopant atom on the Ag-doped Mg surface which shows a very large dissociation 
barrier. 

The dissociation of the H2 molecule is only the first step for the absorption of hydrogen. 
A second fundamental step is the diffusion of the products away from the catalytic site. To 
study this, we performed NEB calculations in which the initial state was the final state of 
the dissociation process, and the final state was obtained by displacing one H into a nearby 
hollow site. Figure [3] shows the diffusion path of one of the two hydrogen atoms on the 
Fe-doped Mg surface as an example. The MEP's for the diffusion processes are also shown 
in Fig. [H We observe that the height of the diffusion barrier ^diff is strongly anti-correlated 
to the height of the dissociation barrier Ediss (see Table UTTI) . In fact, Ti, V, Zr and Ru have 
zero dissociation barriers, but they bind the products very strongly, which results in high 
values of ^diff- By contrast, Ag, Cu and Pd produce large dissociation barriers, but they 
have low diffusion barriers (in fact, no barrier at all for Ag). In between there are Fe, Ni 
and Rh, which represent the best compromise in combining low activation barriers for both 
processes. Ni is the best possible choice overall. 

We note in passing that the catalytic effect of Ni dopant on MgH2 for the dehydrogenation 
process (not studied here) has been experimentally demonstrated by Jensen et al. 64| . 
showing an activation energy reduced by 0.5 eV with respect to that obtained with pure 
MgH2. 

It is interesting to correlate the height of the barriers with the position of the (i-band 
of the transition metal dopant with respect to the Fermi energy Ep (here we define the 
d-band, Pd{E), as the projection of the electronic density of states onto d type spherical 
harmonics). In particular, it is useful to consider the first energy moment of the (i-band, 
or d-band centre, defined as Ed = J^^dE{E — EF)pd{E), where Eq is some cutoff energy 
which we chose to be 7 eV above the Fermi energy. In Fig. H] we plot the dissociation and 
the diffusion energy barriers, Ediss and Edijf, against Ed for all the systems explored (see 



Table [HTj) . It is obvious that the heights of the barriers are strongly correlated with the 
position of the d-band centre. 

The step limiting process in the hydrogen absorption is the one with the largest energy 
barrier between dissociation and diffusion, so the best dopant is the one which minimizes 
the largest energy barrier. It is customary to define the activity of a catalyst in terms of the 
rate of the reaction which is being catalysed. This can often be accurately approximated by 
an Arrhenius relation, and therefore the natural logarithm of the rate is proportional to the 
negative of the activation energy barrier. We can then interpret the maximum of the two 
barriers shown in Fig. [Has indicating the activity of the catalyst. If we draw a line across 
these points, we see that the various transition metals investigated here fit on an inverse 
volcano plot, with Ni, Fe and Rh sitting near the top of the volcano, and therefore being the 
most active catalysts. 

As a matter of interest, in Fig. Owe plot the energy difference between the final and the 
initial state E^^~^^ (both for the dissociation and for the diffusion process) as a function of 
the d-band centre. We can observe some correlation between the two quantities, although 
this is less strong than that observed in Fig. H] for the height of the two energy barriers. If 
follows that the correlation between the energy barriers and E^^~^^ is also weaker than that 
between the energy barriers and the (i-band centre, leaving the latter a better parameter to 
characterize the catalyst. 

From an inspection of Figs. |4]and [5] the d-band centre correlation is evident, and points 
to an ideal d-band centre value of about -1.29 eV. This value cannot be obtained with any 
of the TM-doped Mg surfaces investigated here. Recently, Vegge et al. [9| have investigated 
magnesium 3d TM alloys. They showed that the d-band centre values of the expanded 
alloys obtained with TM belonging to the first raw of the Periodic Table range from +0.93 
to -6.88 eV going from MgSc to MgZn. In particular, MgCu gives a value of -2.37 eV 
while the neighbor MgNi -0.82 eV. It would be interesting to broaden their investigation 
to Ad and TMs to see if the optimal d-band centre value of about -1.3 eV that we have 
extrapolated here could be obtained with some alloys, but this is beyond the purpose of 
the present investigation. 



IV. CONCLUSIONS 



We have performed here a systematic DFT/PBE study of hydrogen dissociation and 
subsequent diffusion over Mg surfaces doped with different transition metals. The dopants 
investigated were Ti, Zr, V, Fe, Ru, Co, Rh, Ni, Pd, Cu and Ag. We have observed that the 
transition metals on the left of the periodic table (Ti, V, Zr), together with Ru, eliminate 
the dissociation barrier altogether, however, the products stick too strongly to the metal 
dopant, therefore hindering diffusion away from the catalytic site. This would result in a 
quick deactivation of the catalyst and therefore a slow absorption process. On the contrary, 
the transition metals on the right of the periodic table do not bind too strongly the H atoms 
(in fact, Ag does not bind them at all), allowing easy diffusion, however, their effect on the 
dissociation barrier is small. We have shown that these two opposite catalytic properties are 



well correlated to the position of the d-band centre, according to the Hammer & N0rskov [65 1 
model. In fact, we have shown that the catalytic activity for the H absorption process can 
be described well by a volcano plot, with the most active catalysts Ni, Fe and Rh sitting 
near the top of the volcano. 
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List of Figures 



1. Minimum energy patlis for tlie dissociation of the H2 molecule, and subsequent diffu- 
sion of one of the two H atoms, on a pure Mg surface and on (Ti, V, Fe, Co, Ni, Cu, 
Zr, Ru, Rh, Pd, Ag)-doped Mg surfaces. 

2. (Colour online) H2 (dark red) dissociation over the Ag-doped Mg surface as viewed 
from side (top figures) and top (bottom figures) positions respectively at IS (left-hand 
panel), TS (central panel) and FS (right-hand panel). The Mg, Ag and H atoms are 
represented respectively by light grey, dark grey and black colours. 

3. (Colour online) H (dark red) diffusion on the Fe-doped Mg surface as viewed from top. 
Figures show positions at the final state of the dissociation which is the initial state 
for the diffusion process (left), at the transition state (centre) and final state (right) of 
the diffusion process. The Mg, Fe and H atoms are represented respectively by light 
grey, dark grey and black colours. 

4. Activation energy barrier for hydrogen dissociation (black) and diffusion (red) of hy- 
drogen on pure Mg and metal-doped Mg surfaces as a function of the rf-band center 
positions. The dashed lines have been drown for eye guidance only. 

5. The energy difference between the final and initial state, E(FS-IS), for hydrogen dis- 
sociation (black) and diffusion (red) on pure Mg and metal-doped Mg surfaces as a 
function of the (i-band center positions. The dashed lines have been drown for eye 
guidance only. 
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